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ABSTRACT 

The effects of late gas accretion episodes and subsequent merger-induced starbursts 
on the plioto-chemical evolution of elliptical galaxies are studied and compared to 
the picture of galaxy formation occurring at high redshift with an unique and intense 
starburst modulated by a very short infall, as suggested by Pipino & Matteucci (2004, 
Paper I). By means of the comparison with the the colour-magnitude relations and 
the [< Mg/Fe >y]-o' relation observed in ellipticals, we conclude that either bursts 
involving a gas mass comparable to the mass already transformed into stars during 
the first episode of star formation and occurring at any redshift, or bursts occurring 
at low redshift (i.e. z<0.2) and with a large range of accreted mass, are ruled out. 
These models fail in matching the above relations even if the initial infall hypothesis 
is relaxed, and the galaxies form either by means of more complicated star formation 
histories or by means of the classical monolithic model. On the other hand, galaxies 
accreting a small amount of gas at high redshift (i.e. z>3) produce a spread in the 
model results, with respect to Paper I best model, which is consistent with the ob- 
servational scatter of the color-magnitude relations, although there is only marginal 
agreement with the [< Mg/Fe >v]-o' relation. Therefore, only small perturbations to 
the standard scenario seem to be allowed. We stress that the strongest constraints to 
galaxy formation mechanisms are represented by the chemical abundances, whereas 
the colours can be reproduced under several different hypotheses. 

Key words: galaxies: ellipticals, chemical abundances, formation and evolution - 
galaxies: merging 



1 INTRODUCTION 

The debate on the formation of spheroids is still one of the 
major issues in astrophysics. In particular, the monolithic 
scenario (e.g. Larson, 1974, Arimoto & Yoshii, 1987, Mat- 
teucci & Tornambe', 1987) and the hierarchical paradigm 
(e.g. White & Rees, 1978, Kauffmann & Chariot 1998, Stein- 
metz & Navarro 2002, Meza et al. 2003) are the main com- 
peting theories on the market. 

According to the monolithic scheme, ellipticals are as- 
sumed to have formed at high redshift as a result of a rapid 
collapse of a gas cloud, characterized by a short and intense 
starburst, followed by a galactic wind at very early epochs. 
This scenario naturally explains the observations of nearby 
ellipticals showing a Mass-Metallicity relation (e.g. CaroUo 
et al., 1993; Gonzalez, 1993; Davies et al., 1993; Trager et 
al. 1998, 2000) and a reddening of the stellar light (Color- 
Magnitude relation, hereafter CMR, Bower et al., 1992) with 
the velocity dispersion of the galaxies, because the most 
meissive ellipticals, having the deepest potential wells, are 
the most metal rich. The lack of significant change in the 



slope and in the scatter of the CMR (Stanford et al., 1998, 
but see Kaviraj et al., 2004), the slow evolution of colors 
(e.g. Sagha et al. 2000; Ellis et al., 1997), mass to light ra- 
tios (e.g. van Dokkum & Franx, 1996; van Dokkum & Ellis 
2003; Rusin et al., 2003; van de Ven et al., 2003) and line 
strength indices of cluster early-type galaxies out to z ~ 1 
(e.g. Bernardi et al., 2003; van Dokkum et al. 1998), are con- 
sistent with a passive aging of the stellar populations and 
hint toward a high and quite synchronized redshift of for- 
mation for these galaxies. Further support to this picture 
seems to emerge from the very recent observations which 
lead Cimatti et al. (2004) to the identification of four old, 
fully assembled, lO^^M© spheroids at 1.6 < z < 1.9. 

On the other hand, hierarchical semi-analytic models 
predict that ellipticals are formed by mergers of already 
formed stellar systems, which closely follow the behaviour of 
their dark matter haloes and regulate the chemical enrich- 
ment of the system (White & Rees, 1978). In this picture 
massive ellipticals form at relatively recent epochs through 
major mergers between spiral galaxies (e.g. Kauffmann & 
White, 1993) occurring over a large redshift interval. Ev- 
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idences favouring the hierarchical scenario are represented 
by the observed interacting galaxies (on-going mergers), el- 
lipticals with disturbed morphologies (i.e. counter-rotating 
cores, dust lanes, ripples, e.g. Kormendy & Djorgovski, 
1989), the morphology-density relation in clusters (Drossier 
et al., 1997), the relatively large values of the Hp index in a 
non negligible fraction of nearby ellipticals (Gonzalez 1993, 
Trager et al 1998). 

We refer the reader to Paper I, Peebles (2002), Calura 
et al. (2004) for a more detailed discussion on the topic and 
references. Here we recall that in Paper I we showed that the 
majority of the optical properties of elliptical galaxies can 
be simultaneously reproduced under the assumption that 
the formation process, occurring at high redshift (as in the 
monolithic collapse scenario), is stronger and faster in more 
massive objects with respect to the less massive ones. In 
fact, the [< Mg/Fe >v] ratio in the cores of ellipticals in- 
creases with galactic mass (Worthcy ct al. 1992; Weiss et 
al. 1995; Kuntschner 2000, Kuntschncr et al., 2001), being 
always larger than zero in the core of bright galaxies (e.g. 
Faber et al., 1992; CaroUo et al., 1993; Davies et al., 1993; 
Worthey et al., 1992). This strongly suggests that the star 
formation (SF) lasted for a period shorter than the time at 
which the Fe pollution from type la SNe becomes dominant 
(see Weiss et al., 1995, Matteucci, 1994). Recently, a simi- 
lar suggestion was proposed by Granato et al. (2004) in the 
scheme of the joint formation of ellipticals and QSOs. 

In order to complete the study, we aim at verifing 
whether Paper I conclusions can be achieved with a more 
complicated formation history, namely to study in a simple 
way the effect of merger induced starbursts on the properties 
of a dominant stellar population assembled at high redshift. 
A complementary view on the analysis made in this paper 
is presented by Thomas (1999) and Thomas et al. (1999). 
In particular, they compared the outcome of a fast collapse 
to the results obtained by the merger of two spirals and 
found difficult to reproduce the [< Mg/Fe >v] overabun- 
dance in the latter case, unless very flat IMFs are invoked 
during the merger. Furthermore, since the [< Mg/Fe >v]- 
(j relation seems to be independent from the environment 
and correlates with galactic age (the older is the galaxy the 
higher is its [< Mg/Fe >v]), Thomas et al (2002) con- 
cluded that SF histories (SFHs) based on the hierarchical 
clustering scenario (Kaulfmann & Chariot, 1998) cannot re- 
produce the observed chemical properties of elliptical galax- 
ies. On the other hand, by means of broad-band photometry 
alone, Shioya & Bekki (1998) constrained the epoch of the 
major merger. They found that mergers between spirals are 
allowed even at z=0.3, but their models cannot account for 
the Mg enhancement. Therefore, it seems necessary to cou- 
ple the information coming from the chemistry to those from 
the broad-band colours in order to better constrain the SFH 
of spheroids. Due to the fact that we are not modelling the 
dynamical evolution, our treatment of the galaxy merging is 
quite simple and focused only to study how the mean chem- 
ical and photometric properties of the stellar populations 
change as a consequence of late episodes of star formation 
triggered by gaseous mergers. In particular, we analysed dif- 
ferent scenarios involving high mass spheroids (i.e. initial gas 
mass Mi^rn = 1O^'^~^^M0), as predicted by Paper I's best 
model, that we consider our fiducial case, and we took it as 
a basis for studying simulated accretion of primordial gas 



triggering a starburst. We explored several cases by varying 
the duration of the starburst and its intensity. 

This paper is organized as follows: in section 2 we 
present the model. In section 3 we discuss the results and in 
section 4 we draw some conclusions. In order to be consis- 
tent with the results of Paper I, we assume a cosmological 
model with Qrn = 0.3, f^A = 0.7 and Ho = 70kms"^ Mpc~^. 
For all models we assume that the star formation starts at 
redshift Zf=5. Variations in Zf do not significantly affect the 
conclusions. 



2 THE MODEL 

The adopted chemical evolution model is based on that pre- 
sented in Paper I. In this particular case, however, we con- 
sider our model galaxies as a single zone extending out to one 
effective radius, with instantaneous mixing of gas. Moreover 
wc take explicitly into account a possible mass flow due to 
the galactic wind and a possible secondary episode of gas ac- 
cretion. Therefore, the equations of chemical evolution take 
the following form: 
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where Gi (t) is the normalized mass density of the element i 
at the time t in the ISM. Xi{t) is defined as the abundance 
by mass of the element i. We address the reader to the orig- 
inal papers of Matteucci & Greggio (1986), Matteucci & 
Tornambe' (1987), Gibson (1997) and Paper I for a compre- 
hensive discussion of these equations. By means of this equa- 
tion we can calculate the evolution of 21 elemental species. 
In particular, the integrals at the right-hand side of the equa- 
tions give the rate at which the element i is restored into the 
interstellar medium both as newly-synthesized and already 
existing element by low- and intermediate-mass stars. Type 
la and Type 11 supernovae (SNe la and SNe II, respectively). 
Here we recall that for SNe la we assume as a progenitor 
model the single degenerate one: a C-0 white dwarf which 
accretes material from a companion (the secondary), even- 
tually reaching the Chandrasckar mass and exploding via 
carbon-deflagration (Whelan & Iben 1973). We make use 
of the SNe la rate formalism derived by Greggio & Renzini 
(1983). The predicted SNe la explosion rate is constrained 
to reproduce the present day observed value in ellipticals 
(Cappellaro et al. 1999). This can be done by means of the 
parameter A in eq. (1), which represents the fraction of bi- 
nary systems in the IMF able to give rise to SNIa. 
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2.1 The initial burst of star formation 

Wc refer to Paper I Model II input parameters for our fidu- 
cial model, namely Salpeter (1955) IMF, Thielemann et al 
(1996) yields for massive stars, Nomoto et al (1997) yields 
for type la SNe and van den Hoek & Groenewegen (1997) 
yields for low- and intcrmcdiatc-mass stars. According to 
this model, the infall timescale (t) and the star formation 
efficiency (u) decreases and increases with increasing galac- 
tic mass, respectively. 



2.2 The merger-induced starburst 

In this paper we need to take into account the possible pres- 
ence of a second infall episode (which simulates the merger), 
therefore we added a merging term in the right-hand side of 
the chemical evolution equation: 

.dGiif). „ , , 

( )acc = Ai,accCe -^""c . {2.) 

Xi^acc describes the chemical composition of the accreted 
gas, assumed to be primordial. This is a good assumption 
for all practical purposes, since, even in the case of gas with 
average Intergalactic Medium metallicity {Z ~ 0.04Zq , e.g. 
Calura and Matteucci, 2004), the metal abundances of the 
accreted gas arc much lower than the ones in the ISM. C 
represents the normalization constant, evaluated by inte- 
grating the law from tacc, namely the time at which the 
merging episode staxts, to the present time, and requiring 
that the galaxy accretes Mace solar masses of gas. Finally, 
Tacc is the gas accretion timescale, fixed to be 0.2 Gyr. The 
reason for this choice is linked to the fact that we expect a 
Tacc ~ T for accreting a gas mass Mace ~ Mium- In fact, we 
recall that r ~ 0.2 — 0.4 Gyr in Paper I's best model. Small 
variations in Tacc, however, do not affect our results. Only a 
much longer timescale {race ~ ^Gyr) would imply that most 
of the newly accreted gas cannot be consumed during the 
second SF episode. 

We recall that, once the conditions for the galactic wind 
onset are satisfied in a given time-stop, we allow the galaxies 
to lose all the residual gas present at that time. Only dur- 
ing the merging-induced starburst, the mass flow rate W{t) 
entering eq. (1) is assumed to be proportional to the star 
formation rate {W{t) = i>{t)) as suggested for starbursting 
galaxies (e.g. Heckman, 2002). It can be shown that these as- 
sumptions do not dramatically alter our conclusions. A more 
detailed treatment for the energetics and the mass exchange 
after the galactic wind has been developed in a separate 
paper (Pipino et al. 2005). 

We limit our analysis to a one-zone treatment out to 
one effective radius, in such a way that the accreted gas is 
uniformly distributed in this region. In this paper, we are 
interested to the effect of star formation histories made by 
multiple episodes, therefore wc force the starburst to happen 
at the beginning of the merger (i.e. at time tacc)- We stress 
that this condition is likely to hold, since we are dealing with 
a large amount of infalling gas, and thus, strong cooling. 

Since we are dealing not only with old stellar popula- 
tions, the light originated by the stellar populations born 
during the second burst could be a significant fraction of 
the total flux coming from the galaxy, thus influencing the 
spectral features. Therefore we calculate the mean stellar 
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Figure 1. Time evolution of abundance ratios in the ISM of a 
model c galaxy. 



abundances by weighting on the V luminosity. In particu- 
lar, we follow Arimoto & Yoshii (1987) and Gibson (1997), 
by means of the photometric code by Jimenez et al. (1998). 
In any case, we confirm the previous findings (Matteucci et 
al. 1998; Gibson 1997; Thomas et al 1999) that for old and 
massive objects the mass weighted averages are very close to 
the luminosity weighted ones, whereas, in the case of more 
recent star formation, the former are systematically ~ 0.2 
dex higher than the latter. 

Once we have defined the parameters governing the 
merging event, it is useful to separate the models analysed in 
different groups, each one representing a different region in 
the space defined by the parameters r , Tacc , Mium , Mace , v- 
A simple classification can be done by considering the 
amount of stars formed during the SF episodes. 

Case 1: Models in which the SFH of the main episode creates 
many more stars than the SF during the merger-induced 
burst. 

Case 2: Models in which the SFH of the main episode is com- 
parable to the SFH during the merger-induced burst. 

Case 3: Models in which the SFH of the main episode is less 
efficient than the SFH during the merger- induced burst. 

We will start from Case 1, showing that those models 
are "perturbations" of the Paper I best model, and in many 
cases predict properties still in agreement with the observa- 
tions (Sec. 3.1). Other models belonging cither to Case 1 or 
2 will be presented in sections 3.2 - 3.4. In particular, we 
verify whether a best model with merger-induced starburst 
which reproduces the properties of ellipticals can be found 
(Sec. 3.2 and 3.3). 
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Table 1. Merger parameters and model predictions 



Model M,„^ tacc [< ^ >vl « ^ >v] l< ^ >vl '^V U-V V-K J-K K. 

(Mq) (Gyr) 



Fiducial 



a (Paper I) 


loii 








0.448 


-0.284 


0.049 


-19.2 


1 


.39 


3,10 


1.03 





b 


lo" 


1 


1 


0.292 


0.288 


0.153 


-19.8 


1 


.45 


3.24 


1.07 


0.94 


c 


lo" 


1 


10 


0.133 


0.946 


0.455 


-20.2 


1 


.13 


2.66 


0.91 


0.54 


d 


lo" 


0.1 


1 


0.391 


0.553 


0.284 


-19.4 


1 


.48 


3.24 


1.07 


0.33 


e 


loii 


0.1 


10 


0.389 


0.330 


0.126 


-19.4 


1 


32 


3.02 


1,01 


0.11 


f (Paper I) 


10l2 








0.433 


-0.129 


0.124 


-22.0 


1 


48 


3,28 


1,09 





g 


10l2 


1 


1 


0.289 


0.372 


0.235 


-22.6 


1 


.57 


3,42 


1,13 


0.26 


h 


10l2 


1 


10 


0.171 


0.843 


0.394 


-23.0 


1 


.23 


2.89 


0.98 


0.74 




10l2 


0.1 


1 


0.415 


0.129 


0.145 


-22.0 


1 


.49 


3.30 


1.09 


0.07 


J 


10l2 


0.1 


10 


0.380 


0.330 


0.176 


-22.2 


1.41 


3.21 


1.07 


0.11 


Closed Box 


k (CB) 


loii 








0.761 


-0.896 


-0.450 


-20.0 





76 


2.73 


0.97 





1 


lo" 


1 


1 


0.524 


0.168 


0.089 


-20.0 


1 


.28 


3.14 


1.04 


0.46 


m 


loii 


1 


10 


0.355 


0.840 


0.335 


-20.5 


1 


.09 


2.69 


0.91 


1.00 


n 


lo" 


0.1 


1 


0.735 


-0.196 


-0.064 


-19.7 


1 


,24 


3.11 


1.04 


0.07 


o 


lo" 


0.1 


10 


0.674 


0.349 


0.031 


-19.7 


1 


12 


2.90 


0,97 


0.08 


P (CB) 


10l2 








0.900 


-0.167 


-0.015 


-22.1 


1 


11 


2,97 


1,00 





q 


ioi2 


1 


1 


0.740 


0.388 


0.143 


-22.5 


1 


.30 


3,26 


1,08 


0.98 




10l2 


1 


10 


0.573 


0.844 


0.279 


-23.0 


1 


.09 


2.77 


0.94 


0.97 


s 


10l2 


0.1 


1 


0.889 


0.080 


0.001 


-22.1 


1 


.13 


3.01 


1.01 


0.09 


t 




0.1 


10 


0.847 


0.344 


0.036 


-22.2 


1 


.11 


2.94 


0.99 


0.14 


Miscellanea 


a+AY87 


loii 








0.614 


0.269 


0.820 


-19.8 


2 


.15 


3.85 


1.24 





b+AY87 


loii 


1 


1 


0.431 


0.855 


1.070 


-20.4 


1 


,41 


3.02 


1.04 


1.04 


c+i/=20 


loii 


1 


10 


0.274 


0.090 


0.221 


-20.9 


1 


.01 


2.69 


0.91 


1.85 




loii 


1 


1 


0.305 


0.081 


0.109 


-19.8 


1.43 


3.21 


1.06 


0.65 




lo" 


1 


1 


0.417 


-0.263 


0.028 


-19.3 


1 


,36 


3.06 


1.01 


0.01 




lo" 


1 


10 


0.293 


0.005 


0.194 


-20.4 


0.99 


2.63 


0.89 


0.75 




lo" 


1 


10 


0.355 


-0.137 


0.088 


-19.8 


1 


,14 


2.85 


0.96 


0.15 


•l+'bursf = "■■'• 


lo" 


0.1 


1 


0.428 


-0.026 


0.098 


-19.3 


1 


34 


3.05 


1,02 


0-01 


<'+tburst = "-^ 


loii 


0.1 


10 


0.426 


-0.216 


0.073 


-19.3 


1 


34 


3.05 


1,02 


0.01 


b+A=0.06 


lo" 


1 


1 


0.372 


-0.012 


0.194 


-20.0 


1 


,51 


3.36 


1,11 


0.36 


u(i/=55) 


ioi2 


1 


1 


0.363 


0.696 


0.388 


-22.2 


1 


,67 


3.48 


1.14 


0.48 


v(,. = 110) 


ioi2 


1 


2 


0.416 


0.839 


0.416 


-21.8 


1 


,48 


3.25 


1.08 


1.67 



Paper I's best model and closed box model without mergers are labelled as Paper I and CB, 
respectively. 



3 RESULTS AND DISCUSSION 
3.1 Perturbing the best model 

We start to explore the parameter space by considering 
merger-induced starburst as perturbations of Paper I's best 
model. In this case the merging-induccd starburst is forced 
to last for tburst = 1 Gyr with the same star formation effi- 
ciency of the initial one (i.e. v = lOGyr"^ and v = 22Gyr~'^ 
for models with initial mass in gas Mium = 10^^ and 
Mium = 10^^ M0, respectively), unless otherwise stated. As 
a first step, in fact, the burst duration and efficiency in con- 
verting gas into stars are chosen in order to have duration 
and intensity similar to those characterizing the main star 
formation episode (i.e. before the galactic wind). In this CEise, 
the SFH during the second burst will be influenced only by 
the gas supply {Mace)- 

We explored very extreme cases, namely we run models 
in which the accreted gas mass is either ~ 10% (i.e. those 
probably leading to Case 1) or ~ 100% (which, according 
to their SFH, will belong either to Case 1 or 2) of the ini- 
tial mass, occurring either 1 or 10 Gyr after the galactic 
evolution has started. Mergers occurring at lower redshifts 
are not allowed, to avoid violation of the constraints on the 
present-day Type II SNe rate in ellipticals; in fact, SNII are 
not observed in local ellipticals (Cappellaro et al. 1999). 

We present our results for the fiducial case in Table 1 
(models a-j), where, each galactic model is identified by a 



letter and the initial luminous gas mass in Cols. 1 and 2, 
respectively. We show in Cols. 3 and 4 the amount of gas 
mass accreted during the merger and the time at which the 
accretion begins, respectively; in Cols. 5-7 we present the 
stellar [< Mg/Fe >v], [< Fe/H >v] and [< Mg/H >v]. 
Finally, in Cols. 8-11 we show the V magnitude and the U-V, 
V-K, J-K colours. 

In the last column we present the ratio R* between the 
mass in stars assembled during the merger-induced burst 
and still alive at the present time, and the stars formed dur- 
ing the initial SF episode. We consider the models in which 
0.5 < < 1.5 belonging to Case 2. Before analysing the re- 
sults, it is important to keep in mind that models accreting 
the same amount of gas [Mace] can undergo different SFHs. 
This is due to the fact that when tacc is low the galaxy 
could be still undergoing a galactic wind. Therefore, the net 
amount of gas available for new SF (i.e. the sum of the gas 
left after the initial burst, the gas restored by stellar mass 
loss and the newly accreted gas minus mass lost for the wind) 
might be smaller than the ISM mass of a galaxy which, af- 
ter the wind phase, passively evolves out to, e.g., tacc ~ 10 
Gyr. In this last case, it is more likely the all the accrctd 
mass Mace can be turned into stars. The second group of 
rows (models k-t) will be discussed in Sec. 3.2. Other model 
belonging either to Case 2 or Case 3 are grouped under mis- 
cellanea and will be discussed in sections 3.3 and 3.4. As 
it can be seen from the entries of this Table, the most evi- 
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Figure 2. Fiducial (full squares) and closed box (empty squares) 
model predictions compared with observations. The mean value 
of the observed colour at a given My is given by the solid line (fit 
to Bower et al.'s, 1992, data). The observed scatter is represented 
by the dashed lines (the off'set from the mean value is taken to be 
±3(7, Bower et al., 1992). Models are labelled according to Table 
1. 



dent effect produced by the merging-induced starburst to the 
fiducial model is a decrease in [< Mg/Fe >v], although the 
ISM abundances increase as a consequence of the increased 
stellar nucleosynthesis. The reason lies in the fact that the 
[Mg/Fe] ratio in the ISM after the galactic wind (and before 
the occurrence of the merger) is dominated by Type la SN 
explosions, which keep its value well below the solar ratio. 
When the merger-induced starburst occurs, most stars form 
out of this gas (thus lowering the total [< Mg/Fe >v]); 
at the same time, the new generations of Type II SNe pro- 
duce a sudden increase in the ISM [a/Fe] ratios, but the 
value reached in this way is lower than the typical abun- 
dance ratio of the ISM during the first ~ 500 Myr (see Fig. 
1). Therefore, the new SNII contribution is not enough to 
make [< Mg/Fe >v] rise again significantly. Since the star 
formation efficiency and the infall time scale for the fidu- 
cial case were tuned to reproduce the color-magnitude, the 
mass-metallicity and the [< Mg/Fe >v']-o" relations (see 
Paper I), the straightforward result is that, at this stage, 
the merging episode worsens the results obtained in Paper 
I. We recall that the star formation rate peak values dur- 
ing either the first or the second burst have similar intensity 
(e.g. ~ 8OA/0j/r~^ in the case of models d) . 
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Figure 3. Fiducial (full squares) and closed box (empty squares) 
model predictions compared with observations. Fit to the data 
(solid) and 3(t boundary (dashed) by Bower et al (1992). Models 
are labelled according to Table 1. 
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Figure 4. Fiducial (full squares) and closed box (empty squares) 
model predictions compared with observations. Fit to the data 
(solid) and Sc boundary (dashed) by Bower et al (1992). Models 
are labelled according to Table 1. 



3.2 Comparison with merging-induced bursts in a 
Closed Box case 

A possible alternative route to reconcile the presence of 
mergers with the observed [< Mg/Fe >v], requires per- 
haps more complex star formation histories. As a first step, 
we run also a closed box model (Table 1, models k-t), with 
the same input parameters of the fiducial case, the only dif- 



ference being that in this case the gas mass is already present 
since the beginning. From the analysis made in Paper I, in 
fact, we know that this case gives very high [< Mg/Fe >v], 
at odds with observations. Therefore we are interested in 
checking if a closed box model plus a second starburst trig- 
gered by a merger (which always lowers [< Mg/Fe >v] in 
the stellar populations), can match the observations. In this 
way, the fiducial and the closed box model allow us to re- 
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Figure 5. Fiducial (full squares) and closed box (empty squares) 
model predictions (entire galaxy) compared with observations 
(within O.lRoff). The dashed lines limit the observed values as 
given by Thomas et al (2002) from data by Gonzalez (1993), 
Mehlert et al. (2000), Beuing et al. (2002). Models are labelled 
according to Table 1. The arrows show that the predictions for 
the whole galaxies become consistent with observations, when the 
radial gradient in [< Mg/Fe >v] (a-s predicted by Paper I's best 
model) is taken into account. 



strict a first region of the parameter space (r varying in the 
interval [0,0.4] Gyr) in which also intermediate cases can 
occur, possibly producing normal ellipticals. 

In order to derive more quantitative constraints we 
compare both fiducial and closed box cases with observa- 
tions in Figs. 2-5. In particular, the fiducial models (either 
a or a 10^^ Mq galaxy) are represented by filled 

squares, whereas the closed box galaxies by empty ones. The 
solid (dashed) lines give the fit (3 a boundary) to the ob- 
served values in the case of the three CMRs, as given by 
Bower et al (1992). In the case of the [< Mg/Fe >v]-Mv 
relation (Fig. 5), the lines bracket the observed region as 
given by Thomas et al (2002). 

The fiducial model still reproduces the CMRs in most 
of the cases when mergers occur; the same happens in the 
semi-analytical models applied to the hierarchical cluster- 
ing scenario, which succeed in reproducing the CMRs (e.g. 
Kauffmann & Chariot, 1998), but fail in reproducing the 
[< Mg/Fe >v]-(j relation, as shown by Thomas & Kauff- 
mann (1999). Care should be taken since our predictions 
are made for a one-zone model, thus neglecting radial gra- 
dients in the colours and in the abundances. Furthermore, 
the observations were made by means of a fixed aperture 
which is in general different from the radius considered here 
(Bower et al, 1992; see also Paper I). In any case, the differ- 
ence between the detailed predictions of a multi-zone photo- 
chemical code (such as those of Paper I) and the ones pre- 
sented here are smaller than the effects produced by either 
the merging-induced starburts or by the assumption of an 
instantaneous initial infall. The J-K colour is not an useful 



tool to discriminate among the models. In general, small (i.e. 
Mace < O.lMium) or high redshift (i.e. tacc < 1 Gyr) pertur- 
bations do not significantly alter the colours, thus predicting 
values well within the observed scatter. Cases c and h are 
ruled out, because they predict U-V and V-K colours which 
are 0.3 mag and 0.5 mag less than the typical value at a given 
Mv, respectively. As expected from the analysis made in 
Paper I, models without infall (i.e. closed box) predict very 
blue galaxies. The agreement with the colours is even worse 
when very recent and massive mergers are allowed (e.g. case 
m). Only one case of massive merger seems to be consistent 
with observations in the whole mass range, namely the one 
in which Mace = Mium and tacc ~ 1 Gyr (models b, g). In 
these particular models, in fact, the merging-induced star- 
burst increases the mean stellar metallicity (and thus makes 
the colours redder) and occurs at a sufficiently high redshift 
that it can be considered a continuation of the main star 
formation episode; in other words, it is like to have a sin- 
gle burst converting ~ 1.5 — 2 x Mium of gas into stars and 
evolving passively since then. 

On the other hand, when analysing the trend of [< 
Mg/Fe >v\ versus My, the presence of a radial gradient 
in the abundance ratio should be considered. In fact, we are 
showing predictions which refer to the whole galaxy, whereas 
the observations are taken within 0.1 effective radius. More- 
over, according to Paper I's best model (see their sec. 3.2), 
a variation of the order of A[< Mg/Fe >v] ~ 0.2 - 0.3 
dex is expected to occur within one effective radius, with 
this ratio increasing outwards. These predictions have been 
recently confirmed by the observed gradients in the galax- 
ies NGC4697 (Mendez et al. 2005) and NGC821 (Proctor et 
al., 2005). However, these results are at variance with the 
mean value of the gradient slope (A[< Mg/Fe >v] ~ 0) ob- 
served by Mehlert et al. (2003) for a large sample of early- 
type galaxies in Coma. This topic is still debated in the 
literature, given its importance in the study of the galaxy 
formation process. At this stage further (either theoretical 
or observational) efforts are needed, which go beyond the 
scope of the present work. In any case, in this paper we as- 
sume the results of Paper I on the radial variation of the 
[< Mg/Fe >v\ ratio. When this happens, both models a 
and /, their predictions in the [< Mg/Fe >v\-Mv plane 
move down in the area where the observations lie (see ar- 
rows in Fig. 5). Under the hypothesis that the gradient is 
unaffected by merging episodes, it is clear from Fig. 5 that 
models with small perturbations are still in agreement with 
the data. In the case of the closed box models, we can see 
that the predictions can be reconciled with observations only 
if a very strong (i.e. A[< Mg/Fe >v\ > 0.4) radial gradient 
in [< Mg/Fe >v\ occurs. However, this is not predicted for 
the CB cases and is not confirmed by observations, which 
show A[< Mg/Fe >v\ < 0.2 (Mehlert et al., 2003, Mendez 
et al., 2005). Among this set of models, only cases m and 
r do not require gradients stronger than those predicted by 
Paper I's best model. However, they are ruled out by the 
photometric analysis. In more realistic cases the models fea- 
turing a merging episode would face another problem that 
might make their disagreement with the observations even 
worse. In particular, we expect that the newly accreted gas 
undergoes stronger cooling in the central region, thus lead- 
ing to longer star formation in the galactic core with re- 
spect to the outskirts. In this case the radial gradient of 
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[< Mg/Fe >v] will be even larger than the one predicted 
by the best model of Paper I and also at odds with observa- 
tional results. Therefore, the predictions for the central zone 
of perturbed models lie well below the observed region, at a 
given magnitude. 

Finally, we remark that all models with a non-negligible 
merger, exhibit very high metallicity ([< Fe/H >] ~ 1, 
[< Mg/H >] ~ 0.4), when the second burst occurs very 
late, at variance with observations (e.g. Kobayashi & Ari- 
moto, 1999, who give [< Fe/H >] ~ —0.3). Because of 
the different ages of the two groups of stellar populations in 
the model galaxies (i.e. the old one, made during the first 
episode of SF, and the young one, made several Gyrs later 
during the second burst), we do not consider safe to obtain 
estimates of the metallicity indices by simply using a cali- 
bration relation, as done in Paper I. Therefore, we avoided 
this part of the analysis. Tantalo & Chiosi (2004, see also 
Maraston & Thomas, 2000) reached conclusions similar to 
ours by analysing the spread in the Hjj — [MgFe] plane and 
the colour evolution. In particular, they found that the pres- 
ence of young populations superimposed to the old bulk of 
stars would be detectable, at odds with observations. On 
the other hand, small merging events producing differences 
in the a-enhancement can explain the observational scatter. 
However, in a very recent paper, Thomas et al. (2004) found 
that a non-negligible fraction of low-mass ellipticals in their 
sample shows evidences of a young (i.e. age < 2 Gyr) and 
metal- rich {[Z/H] ^ 0.2) stellar population which represents 
^ 10 — 20% of their total stellar mass. These values are inde- 
pendent from the environment and the typical [a/Fe] of this 
population are 0.2 dex lower than the bulk of the stars, as 
expected if the young component formed out of gas strongly 
polluted by SNIa (e.g. models c and d). 

3.3 Other star formation histories? 

We complete the analysis of both Cases 1 and 2 by con- 
sidering models in which the SFH is modified by changing 
the second burst duration and its star formation efficiency 
(see Table 1). For example, in case c+tburst~0.5 (i.e. model 
c with tturst=0-5 Gyr), the colours arc still bluer than the 
fiducial CEise and very similar to those predicted by model 
c, although the predictions regarding the chemical evolution 
arc significantly better. In order to improve the agreement 
between the predictions of the fiducial case and the obser- 
vations, we need simultaneously tburst~0.1 Gyr, and either 
a low amount of mass accreted or a low tacc (i.e. high red- 
shift). In these cases, in fact, the time interval during which 
new stars are allowed to form, is short enough to ensure 
negligible perturbations in the final abundance estimates. 
This result holds even in very extreme cases, namely we 
run a model in which Mace = O.lMjum, tburst=0.1 Gyr and 
1/ = lOOOGyr"^, for mergers occurring at tacc < 4 Gyr (i.e. 
z~l). Despite the large value adopted for v, in fact, the 
SFR is only ~ lOOMoyr^^, owing to the small amount of 
gas present in the ISM. This finding is particularly impor- 
tant because, allowing for mass accretion onto the galax- 
ies, but limiting the star formation process to a very short 
period, the typical star formation rate is of the order of 
~ 4O-8OM0yr-^ (- 4-6M0yr-^) for model b+tbur8t=0.1 
{d+tburst=0.1). Therefore, these models cannot convert into 
stars all the gas they accreted, and might end up with a 



cold and massive ISM (even ~ 0.4 — 0.5M;„m in cases with 
Mace = Mium), which is not observed in present-day ellip- 
ticals. Allowing for a very high SF efficiency it is helpful in 
removing all the cold gas and, at the same time, it shows 
the robustness of our conclusions for a large range of SFR. 

Finally, we consider cases in which also the SFH of the 
initial burst deviates significantly from the fiducial model, 
in order to complete the study of the parameter space for 
Cases 2 and 3. We concentrate on cases in which the galaxy 
accretes a gas mass Mace = Mium, since in the previous sec- 
tions we have shown that models with Mace = 0.1 Miy,m can 
be considered as small perturbations (see the ratio R* in Ta- 
ble 1). For what concerns Case 2, the aim is similar to what 
done in Sec. 3.2., namely we try to have a new best-fitting 
model for the [Mg/Fe] vs. mass relation. In particular, we 
explore the cases in which we keep r fixed and change f. Val- 
ues of = 35 — 45 Gyr~^ during both the main SF episode 
and the starburst of model g improve only mildly the agree- 
ment with observations. A significant increase in the pre- 
dicted value for [< Mg/Fe >v] (i.e. 0.363 dex) is attained 
when V = 55Gyr~'^ (model u). In this case, we predict the 
ratio R* to be ~ 0.5, therefore it can be regarded as a model 
similar to g, but with a more intense SFH during the second 
burst. The problem is that, however, the colours start to 
deviate from the observed relations (U-V=1.67, V-K=3.48 
for Mv = —22.2 mag) and the Fe abundance becomes quite 
high ([< Fe/H >v] = 0.700). 

We also tried another exercise, namely we abandoned 
our self-consistent evaluation of the time at which the galac- 
tic wind occurs, based on a detailed treatment of the SN 
remnant evolution, and fixed it a priori. Here we present 
the results for model g, with tg^ — 0.25 Gyr. In this 
case, although the predicted abundance ratios (namely [< 
Mg/Fe >v] = 0.433, [< Fe/H >v] = 0.414) are fairly 
consistent with observations, the colours become too blue. 
In fact, we predict U-V=1.22 and V-K=2.83 for a total V 
magnitude Mv = —21.0 mag. 

This ad hoc model can be seen as a possible way to 
extend model g to Case 2. In fact, owing to the adopted 
self-consistent treatment for the development of the galactic 
wind, model g carmot reach _R, ~ 1. Here, instead, we reduce 
the SF timescale during the initial burst (i.e. tg-w) by a factor 
of ~ 4; therefore we have i?» ~ 1. 

Concerning Case 3, we start by increasing the SF effi- 
ciency during the second burst. For instance see the model 
c+v=20 (i.e. model c with v = 20Gyr~^, which leads to 
Ri. ~ 2) , which improves the agreement between the predic- 
tions regarding the chemical properties, although the colours 
arc still bluer than the fiducial case. Finally, we present 
model V, an illustrative case for a 10^^ M© galaxy in which 
we vary the parameters in order to have a mild and pro- 
longed SF before the onset of the merger-induced burst (oc- 
curring at 2 Gyr). This corresponds to the region of param- 
eter space in which r varies in the interval [0.5, 10] Gyr 
and < 20 Gyr^^ (during the first burst). In particular, we 
chose r = 3 Gyr, v = 2.2Gyr~^ and v = llOGyr"^ during 
the burst. A wind developes only after the starburst. The 
model exhibits [< Mg/Fe >v] = 0.416, U-V=1.48 and V- 
K=3.25 (for a A4v = —21.8 mag), values which arc pretty 
close to Paper I's best model results. Nevertheless, the SF 
produces a very high Fe enrichment in the bulk of the stars 
([< Fe/H >v] = 0.839) which is not observed. In the hy- 
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Figure 6. Star formation rate for some representative cases. Up- 
per panel: Solid: Paper I's best model (/). Dotted: The SFH dur- 
ing the main burst and during the merger have similar inten- 
sity and duration (as in Case 2, model h). Two models in which 
Mace = O.lMjjjm are shown by dashed (model i) and long-daslicd 
(model j) lines, respectively. All these models have the same SFH 
during the initial burst. Lower panel: Solid: The SFH during the 
main burst is more intense than during the merger (model g). 
Dashed: The SFH during the main burst and during the merger 
have similar intensity and duration (as in Case 2, model u). Dot- 
ted: The SFH during the main burst is less intense than during 
the merger (as in Case 3, model v). 



pothesis that this kind of galactic model can be fairly rep- 
resented by a SSP (but see Pipino, Matteucci & Chiappini 
2005, in preparation), this would translate into the mean 
iron Lick index < Fe >= 4.5 (once Worthey's 1994 simple 
stellar populations are adopted, see Paper I). 

The effect of moving the epoch of the mergers to tacc > 



5 Gyr is to worsen the results, since the colours become 
progressively bluer (e.g. V - K =3.13 for My = —22.1 mag). 

Given its R„ ^ 2 and the behaviour of its SFH, to some 
extent, this is a complementary result of what Matteucci 

6 Pipino (2005) have shown for PopIII stars in ellipticals 
(see also Gibson 1996 for a more extended discussion on 
the bimodal SF in ellipticals) by changing the IMF at early 
epochs. In any case, even a short pre-enrichment of the bulk 
of the stellar populations, produce either too red galaxies or 
abundance patterns in disagreement with observations. 

In Fig. 6 we compare the star formation histories of 
some models, chosen as the most representative cases shown 
in the previous sections. In the tipper panel we compare the 
Paper I best model (/) SFH with those from the cases h, i, j. 
In the lower panel, the solid line is taken from model g: the 
main burst clearly dominates over the merger-induced SF 
episode. Case 2 is represented by model u (dashed line): the 
two SF episodes have a very close intensity and shape. They 
are separated by a short episode of galactic wind. Finally, 
a model in which the SF during the merger-induced burst 
exceeds the initial episode (Case 3, model v) is shown by a 
dotted line. 

3.4 Other pcirameters effect 

At this point it is useful to discuss other possible sources of 
degeneracy in our model results. First of all, it should be 
noticed that all the cases but the two Paper I's best mod- 
els (labelled as Paper I in Table 1), fail in reproducing the 
present-day SNIa rate observed in ellipticals (0.18 ± 0.06 
SNu, Cappellaro et al 1999). For example, model c predicts 
an explosion rate of ~ 1 SNu, whereas model b gives 0.48 
SNu. The latter value, however, is only marginally consis- 
tent with the rate of 0.3 ±0.1 SNu given by Cappellaro et 
al. (1999) for ellipticals bluer than the average. The rea- 
son for this disagreement is that the parameter A in eq. 
(1), is set equal to 0.18 in order to have a good fit in the 
case of models a and /. For comparative purposes, we run a 
model in which we use a lower value for A (Table 1, model 
b+A=0.06), namely we chose A=0.06 so that the present- 
day SNIa rate for model b becomes 0.19 SNu. This leads 
to luminosities and colours essentially unchanged, whereas 
the [< Mg/Fe >v] ratio rises toward values more similar 
to Paper I's best model predictions, as a consequence of the 
lower Fe abundance in the ISM when the merger-induced 
starburst occurs. However, a systematic variation of A with 
the properties of the merger, would be in contrast with the 
a-priori definition of A. 

Another possible solution can be a variation either with 
time or metallicity of the fraction of binary systems (e.g. De 
Bonder & Vanbeveren, 2002). It should be remarked, how- 
ever, that, in elliptical galaxies, the metallicity increases very 
fast in the first hundreds of Myr, reaching the solar value at 
the time of the galactic wind, and being higher than solar 
in the wind regime. Therefore, the final (i.e. present-day) 
value for A should be achieved well before the occurrence of 
the second burst of star formation. In any case, we do not 
consider this as a viable solution for reconciling star forma- 
tion histories ruled by random gas accretion episodes with 
the observed features of elliptical galaxies. 

We also tried to change the IMF, by using a typi- 
cal flat one (Arimoto & Yoshii 1987, models a+AY87 and 
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b+AY87). This IMF was ruled out in Paper I beeause it 
predicts a too high metallicity in the stellar populations 
and too high values for [< Mg/Fe >v]- The riicrging- 
induced starburst, when occurring in these models, pushes 
the [< Mg/Fe >v] down to a reasonable value for galaxies 
of that given mass. On the other hand, the Fe and Mg abun- 
dances exceed the solar value by a factor of ten, at variance 
with what is suggested by the observations. The possibility 
for a flat IMF during a merging process was explored by 
Thomas (1999), who concluded that a significantly flat ex- 
ponent (i.e. X < 0.8) is needed to satisfy the constraints on 
[< Mg/Fe >v]. 

Finally, we tested the effect of changing the composition 
of the accreted gas. In particular we focused on a solar com- 
position, as given by Anders & Grevesse (1989). When this 
option is applied to model d, the differences with the predic- 
tions shown in Table 1 are negligible. This is a straightfor- 
ward consequence of the fact that the total mass of the ac- 
creted gas which is turned into stars is very small compared 
to the mass of the dominant stellar population. A stronger 
effect is expected for models in which Mace = Mium- In this 
case, in fact, model c would exhibit [< Mg/Fe >v] = 0.258, 
owing to the great importance of the youngest stellar popu- 
lations in the weighted average. Also the colours become red- 
der (by ~ 0.1 mag), but not enough to match the OMRs. On 
the other hand, the effect of assuming Zacc = Zq on model 
h predictions is milder for the chemical aspect, whereas the 
colours are ~ 0.1 mag redder. This implies that, in any case, 
the time of the occurrence of the merger- induced starburst 
is more important than the chemical composition of the ac- 
creted gas. 



4 CONCLUSIONS 

In this paper we completed the analysis of the photochemical 
properties of elliptical galaxies begun in Paper I, where we 
tested the high redshift formation of ellipticals. We showed 

that a photo-chemical diagnostic, namely the CMRs and the 
[< Mg/Fe >v]-Mv, is needed in order to verify whether 
star formation histories differing from one of the simple 
monolithic one lead to the same good agreement with obser- 
vations. This diagnostic allowed us to quantify Paper I's best 
model response to perturbations, represented by merger- 
induced starbursts. Here we summarise our main conclu- 
sions: 

• Paper I's best model seems to preserve the consistency 

with observations, within the observed scatter, only when 
small merging events occur. These results reinforce Paper 
I findings, namely that the star formation process must be 
more efficient and faster in the more massive galaxies in 
order to reproduce the largest possible set of optical ob- 
servablos. The measured scatter in these relations allows 
only small and early perturbations, i.e. Mace ~ Q.lMiam and 
^acc >3. 

• Either bursts lasting more than 0.5 Gyr and involving a 
gas mass comparable to the mass already transformed into 
stars during the first episode of star formation, or occur- 
ring at low redshift (i.e. z=0.2) are ruled out. Only cases in 
which the burst duration is limited to 0.1 Gyr are marginally 
consistent with observations, if tacc is small. 



• Furthermore, these merging models fail in matching the 
above relations in the case in which the initial infall hypoth- 
esis is relax;ed, and the galaxies form through the classical 
monolithic paradigm (e.g. Larson 1974). Among these mod- 
els, galaxies accreting a small amount of gas at high red- 
shift produce a spread in the model results, with respect to 
Paper I's best model, which is consistent with the observar 
tional scatter of the CMRs, but they fail in reproducing the 
[< Mg/Fe >v]-o- relation. 

• On the other hand, models in which the SFH during 
the merger leads to the creation of many more stars with 
respect to what happened during the initial burst, do not 
match the observational constraints on ellipticals. 

• The possibility of intermediate cases is appealing, but, 
because of the tightness of the observed colour-magnitude, 
mass-metallicity and [< Mg/Fe >v]-a relations, it seems 
to be inconsistent with a systematic presence of merging 
episodes, given the stochastic nature of the accretions along 
the galactic lifetime. 

Finally, we remark that the strongest constraints to the 
galactic formation mechanism and to the occurrence of late 
time merging events are provided by the chemical abun- 
dances. 
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